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Abstract — A  thermal  management  simulation  tool  is 
required  to  rapidly  and  accurately  evaluate  and  miti¬ 
gate  the  adverse  effects  of  increased  heat  loads  in  the 
initial  stages  of  design  in  all-electric  ships.  By  reducing 
the  dimension  of  Navier-Stokes  and  energy  equations, 
we  have  developed  one-dimensional  partial  differen¬ 
tial  equation  models  that  simulate  time-dependent  hy¬ 
drodynamics  and  heat  transport  in  a  piping  network 
system.  Besides  the  steady-state  response,  the  com¬ 
putational  model  enables  us  to  predict  the  transient 
behavior  of  the  cooling  system  when  the  operating 
conditions  are  time-variant.  As  a  demonstration  case, 
we  have  performed  a  thermal  analysis  on  a  realistic 
naval  ship. 

I.  Introduction 

Advances  in  sensors  and  weapons  for  modern  naval 
combatants  are  resulting  in  orders  of  magnitude  increases 
in  power  draw,  which  place  unprecedented  demands  on 
cooling  systems;  this  trend  is  expected  to  continue  for  the 
foreseeable  future.  As  a  result,  the  design  of  the  cooling 
system  can  no  longer  be  left  to  the  later  stages  of  design, 
but  must  be  considered  much  earlier.  A  simulation  tool  is 
required  to  assess  the  impact  of  the  cooling  system  design 
on  the  overall  ship  in  the  early  stages. 

In  this  paper,  we  present  a  system-level  simulation  tool 
that  provides  rapid  design,  modeling  and  analysis  of  chilled 
water  cooling  systems  using  first  principles.  We  begin  with 
a  short  list  of  required  inputs  appropriate  to  the  early 
stages  of  design: 

1)  Basic  ship  geometry  including  ship  length  and 
beam. 

2)  The  locations  of  major  bulkheads  and  decks. 

3)  The  locations,  load  value,  and  vital /non- vital  des¬ 
ignation  of  heat  loads. 

The  software  contains  specifications  of  many  existing 
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cooling  system  components  such  as  chillers  and  heat  ex¬ 
changers,  including  equipment  in  current  use  and  new 
state-of-the-art  components.  The  modular  nature  of  the 
software  also  provides  the  capability  to  insert  additional 
components,  thus  enabling  the  analysis  of  the  impact  of 
newly  developed  technologies. 

Using  the  input  described  above  and  the  user’s  selec¬ 
tion  of  components  and  topology  templates,  the  program 
automatically  provides  a  system  design  which  includes  the 
following  specifications: 

1)  Thermal  zone  extents. 

2)  Piping  path  and  diameter. 

3)  Specifications  and  locations  of  chillers,  heat  ex¬ 
changers,  pumps  and  valves. 

The  framework  of  the  system  modeling  is  based  on  the 
work  of  Sanfiorenzo  [1]  and  Fiedel  [2]. 

Following  system  design,  the  simulation  tool  performs 
comprehensive  static  and  transient  hydrodynamic  and 
thermal  analyses,  and  computes  temperature,  pressure  and 
velocity  across  the  cooling  network.  Our  formulation  is  fast 
to  evaluate  and  is  suitable  for  exploring  various  configura¬ 
tions  in  the  initial  stages  of  design.  As  a  demonstration 
case,  we  have  performed  static  and  transient  analysis  for 
a  naval  ship  with  realistic  specifications. 

II.  Quasi  one-dimensional  transient  model 
A.  Modeling  assumptions 

One-dimensional  models  of  flow  in  elastic  pipes  are 
derived  by  reducing  the  dimensionality  of  the  full  Navier- 
Stokes  equations  under  the  following  assumptions: 

1)  Quasi  one-dimensional  flow:  this  assumption  as¬ 
sumes  that  the  tangential  and  radial  velocity  com¬ 
ponents  are  zero,  i.e.  ur(r,  9 ,  x,  t)  =  ug(r ,  9,  x,  t)  = 
0,  adopting  a  cylindrical  coordinate  system.  As  a 
result  of  this  assumption,  using  the  full  Navier- 
Stokes  equations,  one  can  show  that 

dp  =  dp  = 
dr  89 

2)  Axial  symmetry:  this  assumes  that  any  quantity 
such  as  u(r,  9 ,  x ,  t )  is  independent  of  the  angular 
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location  6,  and  therefore  u  :=  u(r,  x,t).  This 
assumption  becomes  more  accurate  if  the  local 
curvature  of  the  pipe  is  small.  For  instance,  in 
cases  where  flow  undergoes  a  sharp  turn,  such 
as  a  90°  bend  or  a  junction,  the  axial  symmetry 
may  not  be  a  valid  assumption.  In  those  cases,  we 
model  the  component  (i.e.  bend,  junction,  etc.) 
separately. 

3)  Fixed  radial  dependence:  the  axial  velocity  profile 
is  assumed  to  be  in  the  form  of 

ux(r,x,t)  =  u{x,t)g(r) 

where  u(x,  t )  is  the  mean  of  the  profile  and  thus 
/o  '  Trg(r)rdr  =  irr2,  where  ri  is  the  internal  radius 
of  the  pipe.  The  choices  for  radial  dependence 
functions  depend  on  the  flow  condition.  For  in¬ 
stance,  the  Poiseuille  flow  profile  may  be  assumed 
for  a  laminar  flow  condition,  or  an  experimental 
profile  may  be  considered  for  a  turbulent  flow 
condition.  As  will  be  demonstrated  later  in  the 
paper,  the  explicit  profile  of  g(r)  does  not  appear 
in  the  mathematical  modeling.  However  g(r )  di¬ 
rectly  affects  the  friction  factor  as  different  choices 
of  g(r)  determine  the  shear  stress  at  the  pipe  wall. 

4)  Incompressible  flow:  we  assume  that  the  working 
fluid  is  incompressible.  Therefore,  the  internal 
energy  e  is  expressed  as: 

e  cpT, 

where  cp  is  the  specific  heat  capacity. 

5)  Negligible  kinetic  energy:  the  kinetic  energy  is 
neglected  as  it  has  significantly  smaller  values 
than  the  internal  energy  of  liquids  such  as  water. 
Thus  e  +  u2  / 2  ~  e. 

6)  Negligible  conduction:  in  most  realistic  cases,  ax¬ 
ial  and  radial  conduction  terms  are  negligible 
compared  to  the  advection  terms. 

The  first  three  assumptions  mentioned  above  are  anal¬ 
ogously  used  for  temperature,  and  thus  temperature  can 
be  expressed  as  T  :=  T(x ,  t)h(r ),  where  T(x,  t)  is  the  mean 
temperature. 


B.  Mathematical  Modeling 


1 )  Hydrodynamics:  To  model  the  propagation  of  tran¬ 
sient  flow  in  pipes,  both  the  compressibility  of  the  liquid 
and  the  elasticity  of  the  pipe  must  be  modeled.  Local 
increase  in  fluid  pressure  results  in  local  enlargement  of  the 
pipe  cross  section  area.  In  general  the  unknown  variables 
are  pressure  p(x,t),  axial  velocity  u(x,t )  and  the  cross 
sectional  area  A(x,t).  A  purely  elastic  model  provides  a 
pressure-area  relation.  The  Laplace’s  law  provides  such  a 
relation:  __  _ 

P  =  Pext  +  ~  77),  (1) 


where 


P  ~ 


2  pc2 

71 


with 


K 

!  +  (§)(?)' 


E condor  —  0 


E cond,r  —  0 


Fig.  1  :  Control  volume  for  coolant  water  for  an 
element  inside  the  pipe. 


In  the  above  equation,  c  is  the  speed  of  an  acoustical  wave 
through  the  pipe,  e  is  the  wall  thickness  of  the  pipe,  K  is 
the  bulk  modulus  of  elasticity  of  the  fluid,  and  E  is  the 
Young’s  modulus  of  the  elasticity  for  the  wall  material.  To 
close  the  system,  two  other  equations  are  required,  which 
are  provided  by  the  conservation  of  mass  and  momentum: 


dA  duA 
dt  dx 
du  du2 / 2 
dt  dx 


1  dp 
p  dx 


fu2 

2D 


—  g  sin  9 


(2) 

(3) 


In  the  above  equations,  /  is  the  friction  factor  in  the 
pipe  and  depends  on  the  state  of  the  flow  (i.e.  laminar  or 
turbulent)  and  the  roughness  of  the  pipe  wall.  The  variable 
9  is  the  angle  between  the  pipe  axis  and  the  horizon. 


2)  Energy:  Following  our  assumptions,  we  consider  the 
control  volume  as  shown  in  Figure  1.  The  control  vol¬ 
ume  expands  to  the  pipe  wall  boundary,  where  energy  is 
exchanged  through  conduction.  Moreover,  energy  crosses 
the  boundary  of  the  control  volume  by  means  of  axial 
conduction  and  advection.  The  balance  of  energy  for  a 
generic  control  volume  can  be  written  as: 

£  pedV  =  j>  EadvdS  +  ji  EconddS  (4) 

As  we  mentioned  in  the  previous  section,  the  axial  and 
radial  conductions  are  negligible.  Thus:  §s  EconddS  ~  0. 

The  axial  advection  for  the  control  volume  shown  in 
Figure  1  is  given  by: 


EadvdS  —  Eadv,  x  Eadvx-\-dx 


(5) 


=  pu(e  +  u2 /2)AX 

-  |pu(e  +  u2 /2)AX  +  7-  [pu(e  +  u2 /2)Ax]dx  j> 
d 

=  ~  da:  7l(e  +  “V2)^]  dx 
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Therefore  the  energy  equation  becomes: 


where  T  =  TA. 


fdr 

pc^ai + 


d(uY) 

dx 


)=0. 


(6) 


[pipe#  7|  [At,  uu  rl]/cQmiOl_F  2,  U2,  r2|  [pipe#  2| 

U  >  0 


The  system  of  partial  differential  equations  for  hydro¬ 
dynamics  and  heat  transfer  can  be  expressed  as: 


where: 


dU  H0U 

dt  dx 

F 

(7) 

(A\ 

U=(“J- 

(8) 

/  u  A 

°\ 

H=  huk  « 

0 

(9) 

and 


T  ut 


F=  ini 

V  o 


(10) 


Given  that  A  >  0,  the  matrix  H  has  three  real  eigen¬ 
values.  Thus  the  above  system  can  be  diagonalized  to  the 
following  form: 

aW+A?  =  S  (11) 


dt 


dx 


where: 


(W\ 

W  =  \W2 

\w3 

with  W\ ,  W2  and  W3  being  the  Riemann  invariants  given 
by: 

Wi(H)  =  u  +  4(c  -  c0)  =  u  +  4«/ ^(A1/4  -  Aj/4)  (12) 

y  ^ P 

W2(H)  =  u  -  4(c  -  Co)  =  u  -  4  J  ^(A1/4  -  Aj/4)  (13) 

y  ^ P 

W3(H)  =  Y/A  (14) 

and  A  is  a  diagonal  matrix  whose  entries  are: 

Ai(H)  =  u  +  c;  A2(H )  =  u-c-  A3(H)  =  u.  (15) 


In  most  piping  network  systems,  c  >>  «,  and  therefore 
u+c  >  0  and  u  —  c  <  0.  As  a  result  the  characteristic  curve 
that  is  obtained  from  dx/dt  =  u  +  c  is  forward  traveling 
and  the  characteristic  curve  with  dx/dt  =  u—c  is  backward 
traveling.  The  characteristic  curve  obtained  by  dx/dt  = 
u  can  be  either  forward  traveling  or  backward  traveling 
depending  on  the  sign  of  u. 


C.  Component  Modeling 

1)  One-to-One  Connector:  A  One-to-One  Connector  is 
a  component  that  links  two  pipes  and  provides  boundary 
conditions  for  the  end  of  the  connecting  pipes.  Pumps, 
valves,  bends  and  heat  exchangers  are  such  components. 
These  components  are  modeled  as  zero-dimensional.  The 
schematic  of  the  component  and  the  connecting  pipes  is 


Fig.  2  :  The  schematic  of  the  one-to-one  compo¬ 
nent  model. 


shown  in  Figure  2.  The  six  unknowns  are  shown  in  red. 
These  variables  are  A3,  u\  and  Ti  from  pipe  #1,  and  A2, 
w2  and  r2  from  pipe  f/2.  The  hydrodynamic  variables  Ai, 
A‘2 ,  Hi,  and  u2,  are  independent  of  temperature  and  they 
can  be  solved  separately.  Thus  we  need  four  equations  to 
close  the  system.  Forward  traveling  characteristic  waves 
carry  information  to  the  component  through  pipe  //  1,  and 
backward  traveling  characteristic  waves  carry  information 
to  the  component  through  pipe  ff  2.  These  two  Riemann 
invariants  provide  two  independent  equations.  The  con¬ 
servation  of  mass  provides  the  third  equation.  The  fourth 
equation  is  obtained  by  incorporating  the  total  pressure 
change,  denoted  by  A Pt,  for  fluid  traveling  from  pipe  #1 
to  pipe  f/2.  In  summary  these  four  equations  are: 


=  Wi 

(16) 

=  W3 

(17) 

A\U\  =  A2U2 

(18) 

(P2  +  -ul)  -  (Pi  +  -«?) 

=  APt 

(19) 

The  pressure  change  depends  on  the  component.  For  in¬ 
stance,  for  pumps,  A Pt  represents  the  pressure  increase 
imparted  by  the  pump,  and  thus  A Pt  :=  Pp(Q)  where 
Pp{Q)  is  the  pump  characteristic  curve  as  a  function  of 
volumetric  flow  rate  Q  =  A3ui  =  A2«2.  For  valves,  bends 
and  heat  exchangers,  APt  represents  the  pressure  loss 
as  flow  travels  through  these  components  and  thus  it  is 

generically  expressed  as:  A Pt  =  —  —  pu\u\,  where  K  >  0 
is  the  loss  coefficient. 


Note  that  the  pressure  is  algebraically  related  to  area 
through  the  constitutive  model  given  by  Equation  1. 


To  solve  for  T 1  and  T2  we  use  the  balance  of  energy  and 
the  characteristic  equations.  The  balance  of  energy  for  the 
component  can  be  written  as: 

Qt  =  rncp{T2  -  T3) 

In  this  setting,  Qt  is  the  amount  of  total  heat  transfered  to 
the  component.  Substituting  for  m  =  pu\A\  =  pu2A2  and 
the  definitions  of  Ti  =  TfAi  and  T2  =  T2A2,  we  arrive  at: 

Qt  =  pcp(u2Y2  -  uiTi) 

To  close  the  system,  we  use  the  information  that  propa- 
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gates  along  the  characteristic  dx/dt  =  u.  This  follows: 

if  «>0,  7j  =  W3 
if  u  <  0,  T2  =  W3. 

We  assume  that  pumps,  valves  and  bends  are  adiabatic 
and  therefore  Qt  =  0  for  these  components.  For  heat  ex¬ 
changers,  i.e.  heat  loads  and  chillers,  Qt  is  either  calculated 
from  a  reduced  order  model  or  it  is  assumed  to  be  given. 

2)  One-to-Two  Connector:  One-to-two  connector  links 
three  pipes  at  a  junction.  The  modeling  of  a  junction  is 
similar  to  the  one-to-one  connector  with  nine  unknowns. 
The  junctions  could  be  of  splitting  or  merging  type.  For 
more  details  on  the  modeling  of  the  junctions  see  reference 
[3], 

D.  Numerical  Method 

We  use  Discontinuous  Galerkin  (DG)  method  to  dis¬ 
cretize  the  system  of  partial  differential  equations  given 
by  equation  7.  We  split  each  pipe  to  Ne  elements  and 
within  each  element  Legendre  polynomials  are  used  as  the 
basis.  We  use  upwinded  flux  that  propagates  information 
between  the  elemental  regions  and  the  bifurcations  of  the 
system.  The  equations  for  components  form  a  nonlinear 
system  of  equations  that  are  solved  at  every  time  step 
using  the  Newton-Raphson  method.  An  Adams-Bashfortli 
scheme  is  used  for  the  time  integration.  A  brief  outline  of 
the  numerical  method  introduced  in  [3]  is  presented  here. 
For  more  details  on  the  numerical  method  see  [4]. 

III.  Demonstration  case 


TABLE  I:  Lumped  values  of  heat  loads  and  their 
locations. 


Identifier 

Q{kW) 

x{m) 

y(m) 

z(m) 

1 

532 

-66.67 

-0.34 

6.76 

2 

5050 

-40.00 

0.00 

16.57 

3 

4723 

-13.34 

-0.08 

18.27 

4 

2611 

13.34 

0.43 

16.95 

5 

153 

40.00 

-2.00 

5.41 

6 

378 

66.67 

2.00 

0.01 

B.  Design  variables 

1 )  Thermal  zones:  Six  thermal  zones  are  considered, 
and  two  chillers  are  placed  in  each  zone.  The  zonal  bound¬ 
aries  and  chillers  are  shown  in  Figure  3. 

2)  Pipe  configuration:  A  simple  loop  layout  is  consid¬ 
ered  for  the  supply  and  return  header;  see  [1]  for  more 
details  and  other  options.  The  port  and  starboard  main 
piping  heights  are  5.20  (m)  and  10.20  (nr),  respectively. 
The  offset  between  the  supply  and  return  header  is  0.50 
(m). 

The  supply  and  return  branch  piping  connect  each  heat 
load  to  the  corresponding  headers.  We  model  each  branch 
as  a  two-pipe  segment  with  a  single  90-degree  bend  in  the 
(y  —  z)  plane;  all  flow  in  the  x  direction  occurs  in  the 
headers. 


In  this  section,  we  use  the  computational  tool  to  an¬ 
alyze  the  performance  of  a  cooling  pipe  network  aboard 
a  modern  naval  ship.  In  the  following  we  describe  the 
problem  inputs  and  the  design  parameters  and  then  we 
discuss  the  results. 

A.  Inputs 

1 )  Geometric  parameters:  The  overall  length  of  the 
ship  is  Lg  =  167.29 (to),  and  the  overall  beam  is  Lg  = 
20.34 (to).  We  place  the  coordinate  system  at  the  center 
of  the  ship  with  the  longitudinal  coordinate,  x ,  expanding 
from  —Ls/2  to  Lg/ 2. 

2)  Heat  loads:  There  are  26  heat  loads  aboard  the  ship. 
We  have  considered  the  maximum  values  that  the  heat 
loads  can  acquire  in  different  modes  of  operation.  For  the 
sake  of  simplicity,  we  have  divided  the  ship  into  six  equal 
longitudinal  segments,  or  zones,  and  lumped  all  the  loads 
that  lie  in  each  zone  into  a  single  load.  The  value  of  the 
lumped  load  is  equal  to  the  summation  of  all  heat  loads 
that  lie  within  that  segment,  as  shown  in  Table  I. 

3)  Physical  parameters:  The  Young’s  modulus  of  the 
pipe  is  E  =  1.84  x  1011(F’a),  and  the  ratio  of  the  pipe 
diameter  to  its  thickness  is  assumed  to  be  D/e  =  8  for 
all  pipes.  The  bulk  modulus  of  elasticity  for  water  is  K  = 
2.39  x  1011  {Pa).  The  density  of  water  is  p  =  1028 (kg/m3). 
The  flow  regime  in  all  pipes  is  assumed  to  be  turbulent 
and  the  friction  factor  is  set  to  be  constant  at  /  =  0.04  for 
all  pipes. 


3)  Pipe/pump  sizing:  The  pipe  diameter  for  each 
branch  is  proportional  to  the  corresponding  heat  load,  Q. 
As  an  initial  estimate,  we  use  the  parametric  equation 
introduced  in  reference  [2]: 


D  =  { 


4  KQ 
Cn  ' 


,0.4 


(20) 


The  diameter  is  then  rounded  to  the  closest  standard  sizing 
for  pipes.  The  diameter  of  the  main  supply/return  header 
is  Dh  =  14  (inches).  The  supply  pumps  are  located  at 
the  exit  of  each  chiller.  Therefore  we  have  a  total  of  twelve 
pumps.  The  pumps  are  all  identical  with  the  characteristic 
curve  obtained  from  curve-fitting  the  empirical  data  as 
shown  in  Figure  4. 


C.  Results 

We  assume  that  the  chilled  water  temperature  is  always 
kept  constant  at  Tchm  =  279.817A'  =  44 F.  This  can  be 
achieved  by  setting  the  temperature  at  chiller  outlet  at 
279.817A' .  Since  we  assume  that  the  pipes  are  adiabatic, 
we  expect  to  preserve  this  temperature  in  the  supply 
pipes.  In  this  fashion,  we  can  measure  the  cooling  capacity 
required  at  each  chiller  to  always  maintain  the  supply 
chilled  water  temperature  at  T  =  279.817 K.  We  initialize 
the  temperature  of  all  pipes  at  t  =  0  to  T  =  279.817 K. 

In  Figure  5,  the  coolant  mass  flow  rate  of  three  chillers 
versus  time  is  shown.  Chillers  CH\  and  CH 2  both  belong  to 
the  same  thermal  zone,  but  due  to  the  asymmetry  of  the 
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Fig.  3  :  Double  main  piping  system  with  branches  (plan  view).  Zonal  boundaries  are  shown  with  dashed  green 
line,  the  supply /return  to  heat  loads  are  shown  with  gray  line. 


Fig.  4:  Pump  characteristic  curve  chosen  for  the 
coolant  network.  The  circles  show  the  empirical 
data  provided  by  the  pump  manufacturer.  The 
solid  line  shows  a  parabolic  curve  fit  to  the  data 
points. 


Fig.  5  :  Chilled  water  mass  flow  rate  at  three 
different  chillers  versus  time. 


pipe  network,  these  two  chillers  have  different  mass  flow 
rates.  All  three  chillers  show  convergence  to  a  steady  state 
after  a  transient  period.  The  “high-frequency”  oscillations 
observed  in  the  mass  flow  rate  signal,  are  the  result  of 
acoustic  wave  propagation  throughout  the  pipe  network. 
Since  c  »  u,  these  oscillations  appear  on  much  smaller 
time  scales  than  the  inertial  time  scale  of  the  flow,  i.e. 
D/u. 

The  primary  objective  of  the  current  simulation  is  to 
compute  how  much  cooling  capacity  each  chiller  requires 
to  remove  the  heat  from  the  loads  in  the  steady  operation. 
In  Figure  6  the  instantaneous  cooling  capacity  of  all  twelve 
chillers  are  shown.  For  the  sake  of  clarity,  the  plot  is  divided 
into  three  plots  with  four  chillers  shown  in  each  plot.  The 
cooling  capacity  shown  in  Figure  6  represents  the  amount 
of  heat  that  must  be  removed  by  each  chiller  to  maintain 
the  supply  water  temperature  at  Tcjnu .  As  it  is  shown  in 
Figure  6,  each  chiller  experiences  a  sudden  rise  in  Q ,  after 
some  initial  interval,  which  is  due  to  the  time  it  takes  for 
the  return  water  to  travel  from  the  heat  load  back  to  the 
chillers.  However,  after  the  sudden  rise,  the  demand  for 
cooling  capacity  settles  to  the  steady-state  values. 

Comparing  the  the  total  cooling  capacity  of  all  chillers 
versus  the  total  heat  load  provides  a  robust  verification 
check  for  the  numerical  computation  carried  out  in  this 
section.  In  Figure  7,  the  total  cooling  capacity  of  the 
cooling  network  is  compared  with  the  total  heat  load.  The 
total  heat  load  is  obtained  by  adding  up  the  heat  loads 
given  in  table  I.  This  value  is  Qioad  =  13.35(MW).  It  is 
clear  that  the  total  cooling  capacity  converges  to  this  value. 


D.  Analysis 

The  highest  capacity  is  Q  =  1.66  {MW)  and  it  belongs 
to  the  third  chiller.  It  is  often  preferred  to  choose  identical 
chillers  for  all  thermal  zones  due  to  lower  installation 
and  maintenance  costs.  Thus  for  the  current  ship  we 
would  need  twelve  chillers  with  the  cooling  capacity  close 
to  1.66(MW),  which  is  an  order  of  magnitude  increase 
compared  to  conventional  ships.  A  state-of-the-art  chiller 
with  the  cooling  capacity  of  1.5 (MW),  has  a  volume  of 
5.00(m)  x  2.40(m)  x  2.30 (to)  and  it  weighs  20000  kg.  (See 
for  instance  York,  Sea  Water  Cooled  Centrifugal  AC). 
Therefore  the  total  volume  and  weight  of  the  chillers 
roughly  become  331  m3  and  24  tons.  Moreover,  a  ship  with 
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Fig.  6:  Cooling  capacity  for  all  twelve  chillers  ver¬ 
sus  time.  For  clarity,  only  four  chillers  are  shown 
in  each  plot.  The  coolant  capacity  is  required  to 
maintain  the  chilled  water  at  T  =  Tchm-  The 
cooling  capacity  eventually  reaches  steady-state 
values.  The  steady-state  capacities  are  different  for 
different  chillers,  and  it  varies  from  0.88  MW  at 
CH12  to  1.66  MW  at  CH3. 


Fig.  7:  The  total  cooling  capacity  of  all  twelve 
chillers  compared  with  the  total  heat  load  of  Q  = 
13.45(MW). 


the  current  design  would  suffer  from  high  vulnerability  as 
CH-i  and  CHg,  and  CUg  and  CH io  are  unable  to  sustain  the 
heat  generated  within  the  thermal  zones  they  are  located 
in.  The  current  design  can  only  handle  all  the  heat  loads 
with  all  chillers  operating  concurrently. 

IV.  Conclusion 

We  developed  a  computational  tool  that  is  suitable  for 
the  early  stages  of  design.  Applying  our  computational  tool 
to  evaluate  the  performance  of  a  proposed  cooling  network, 
we  highlighted  some  of  the  challenges  in  using  the  state- 
of-the-art  technology  for  cooling  modern  naval  ships  with 
high  heat  loads. 
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